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SUMMARY 

This report describes a three-dimensional object simulation 
program which has been developed. Input to the program is a 
description of the object in terms of polygons, and light 
source-object-observer geometry. Output from the program is a 
two dimensional array of numbers representing the brightness 
values of the object as would be seen by the observer. Shadows 
are included in the output image. Examples are shown. The 
procedure for simulation of an optical system using the simulated 
images as input is also given, along with examples. 
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OBJECT SIMULATIONS 


Introduction 

A computer program has been developed at the Visibility 
Laboratory for simulating three-dimensional objects. The program 
produces a two dimensional array of data values which represents 
the object's relative brightness or radiance values for a given 
light source-object-observer viewing geometry. An important and 
unique feature of the program is that shadows are realistically 
included in the image. This program has been used to produce 
images for system simulation and mensuration studies. This 
report briefly describes the method used in the program and gives 
examples of both object simulation and system simulation. 

Description of Program 

The ability to computer simulate three-dimensional objects 
12 

is not new ’ . The distinguishing feature of the program being 
presented here, called NMMV, is its ability to realistically 
handle shadows for an arbitrarily placed light source, and its 
ability to separate the hidden surface algorithm from the 
picture-making algorithm. 

In the previous simulation programs used at the Visibility 
Laboratory, the light source was assumed to be coincident with 
the observer. Further, one had to select the picture parameters, 
such as resolution and region of the object to be displayed in 
advance, since the Watkins hidden surface algorithm Involves a 







scan line-by-scan line approach. The NMMV method implements the 
hidden surface algorithm first and uses a separate algorithm to 
generate the pictures afterward. 

The NMMV method is in principle very simple. Objects are 
defined in a local 3-D coordinate system. The object must be 
broken into planar polygons of 16 sides or less. For curved 
objects, the finer the piecing the better the resulting 
approximation. The polygons are defined by specifying the 
position of their vertices (referred henceforth as nodes) in 
3-space, along with specifying a clockwise connection pattern, 
which will be made clear by a simple example. 

See Fig. la. The polygon shown is actually a square. If we 
are to stand upon the first node and walk along the outside of 
the polygon in a clockwise fashion (from the observers point of 
view) we would encounter the nodes of the polygon in the order 1 , 
4> 3 , 2. Thus we would define the polygon as the connection of 
nodes 1, 2, 3» and 4 in the order of 1, 4,3, 2. 

Definition of the polygons with proper clockwise orientation 
is mandatory. As a rule, if one wishes a polygon to be visible, 
it must be seen to be a clockwise polygon. Note that an observer 
from the opposite position in the example above would see the 
polygon as counterclockwise, and hence for that observer the 
polygon would appear invisible, by definition. One wishing to 
define the very thin rectangular panel on Fig. 1b who wants the 
panel to be visible from all angles would have to define both 
sides of the panel to be two distinct polygons, e.g., polygon 
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A = 1, 4, 3, 2 and polygon B = 1, 2, 3, 4. The program 

interprets polygons which are oriented counterclockwise to be 
pointing away from the observer. These polygons are treated as 
if they were invisible to the observer because it is assumed that 
in simulating solid objects the observer always sees the face of 
a panel and never the back. 

The basic algorithm for taking a polygon-approximated image 
and determining its 2-D projection with shading is as follows: 

(1) The polygon approximation is converted to 
the coordinate system of the light source. 

This is accomplished by having the user 
input azimuth and zenith position angles of 
the object and light source with respect to 
an arbitrary coordinate system centered at 
the observer. The program then performs the 
mathematical calculations required to 
express each node of the object in the 
coordinate space of the light source. 

(2) Taking each polygon in turn, the program 
decides if the polygon is visible, i^£., 
still clockwise oriented. If so, the 
program puts it in a visible list and will 
process it further. Otherwise, it places 
the polygon in a hole list, and no further 
processing will be done on it until step 
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(3). Of the polygons in the visible list, 
some will be in front of others and either 
partially or completely prevent light from 
reaching the more distant polygons. An 
algorithm is employed which for any two 
polygons will decide which polygon is in 
front (if either) and which casts a shadow 
on the other (if such is the case). The 

basic method simply involves projecting both 
polygons on a common plane and testing all 

segments of one against the other to 

determine if there are any non-trivial 
intersections. 

If the two polygons do not intersect, both 
remain in the visible list and testing is 
carried out on the next pair of polygons 
(which pair may involve one of the polygons 
from the previous pair). If an intersection 
is found, the polygon upon which the shadow 
is cast is divided into visible and shadowed 
components. That is, the polygon is 

subdivided into two or more sub-polygons, 
each of which have the property of being 
either completely visible (as far as the 
shadowing polygon is concerned) or 
completely blocked. 

T 
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The sub-polygons which are shadowed are 
transferred to the hole-list and will not be 
processed further until step three. The 
visible sub-polygons are returned to the 
visible list and are thenceforth treated as 
any other visible polygon in that list. 
Thus, they may be tested again against other 
polygons, though they will not be tested 
against their siblings nor will they be 
tested against any polygons their parent 
polygon was already tested against. 

In this manner, all visible polygons are 
tested against each other until the visible 
polygon list is purified and there remains 
in the list only polygons which are 
completely visible to the light source. 
Conversely, those polygons in the hole list 
are either completely shadowed or else were 
oriented in the wrong direction and would 
have been shadowed if they had been put to 
test anyway. This completes step two. 

(3) Now the polygons are all put together in a 
new list, reforming the object more or less. 
As each polygon is put in the new list it is 
tagged as having been either visible or 
not-visible (shadowed). The object is then 
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redefined in the coordinate system of the 
observer, by suitable coordinate 

transformation similar to those used in step 
( 1 ) . 

(4) Polygons are now tested to see if they are 
clockwise. If they are not, they are 
pointing the wrong way and thus must be 
blocked by some polygon(s) pointing the 
correct way. Thus, the observer cannot 
possibly see them and we put them in a new 
hole list which is the same as throwing them 
away since we will never be doing any more 
processing on them. Those polygons which 
have a chance of being visible are put in 
the visible list and will be processed 
further in step (4). It is important to 
remember that of these visible polygons, 
some will have been tagged as lighted and 
some as unlighted since some may have been 
invisible to the light source but are 
visible to the observer. Some will have 
been visible to the light source but will 
not be seen by the observer, and some will 
be visible to both. 

(5) It is the purpose of step five to carry out 
this final partitioning. The method is 
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exactly analogous to that used in step ( 2 ) 
above, so only the results will be discussed 
here. That is, step (5) results in a list 
of pure visible polygons in the sense that 
all polygons remaining in the visible list 
after step four are completely visible to 
the observer. Further, all polygons which 
are in the hole list are blocked from the 
view of the observer. The important 
difference between steps (5) and (2) is that 
the visible and lighted list and a visible 
but shadowed list such that all members of 
the former list are completely visible and 
completely lighted, and all members of the 
latter list are completely visible but are 
not illuminated by the source. 

Finally, the polygons are shaded according 
to their reflectivities and angles between 
the light source, object and observer. That 
is, each polygon is assigned an intensity 
value between one and zero. Shadowed 
polygons may receive non-zero values if the 
user has selected an option of diffuse 
lighting. This allows soft shadows to be 
produced. 
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(7) A separate program can then be used which 
fills in an array of arbitrary size by 
assigning each pixel in the array the 
intensity of the polygon (if such exists) 
which occupies that position in the 

projected plane perpendicular to the viewing 
angle of the observer. Thus, once the basic 
algorithm has been run, the user may choose 
any size for display, whereas in previous 
programs, the size had to be fixed in the 
beginning. 

Examples of Object Simulation 

Figure 2 shows examples of a simulated object. The light 
source is behind the observer and below the horizon. The object 
is shown at azimuth angles of -20, -10, 0, 10, and 20 degrees 

with respect to a plane through the observer and the sun. Zenith 
angle of the object is approximately 30°. The object is rotating 
approximately 20° from one position to the next in a direction 

the same as the apparent rotation due to azimuth angle. Thus the 
total apparent rotation to the observer is about 30° from one 

position to the next. The object is shown with both soft and 

hard shadows. Hard shadows are those areas that receive no 
direct illumination from the source. Soft shadows are the hard 
shadow areas filled in with low brightness values, still 
proportional to the reflectance of the panels in the shadowed 
areas. This serves to show the true geometrical shape of the 
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object. It could also simulate in a crude way illumination of 
shadowed areas by reflected light within the object or from 
another object. 

System Simulation 

The object modeling program provides the input to a system 
simulator which then applies the characteristics of a given 
sensor system. One such system which has been simulated is a 
system employing adaptive optics for atmospheric turbulence 
compensation. (Ref. 3) The following steps are involved in the 
simulation of a long time-exposure image. 

(1) Fast Fourier transform the object radiance 
map. 

(2) Compute the telescope diffraction 
modulation transfer function. 

(3) Compute the modulation transfer function 
associated with the residual wavefront 
errors of the adaptive optics system. 

(4) Compute the modulation transfer function 
associated with miscellaneous errors. 

(5) Multiply the object Fourier transform by 
all these transfer functions, and inverse 
transform to get the degraded image, 
without sensor effects. 





(6) Scale the degraded image spatially and in 

intensity so that the value at each pixel 
corresponds to the mean number of 

photoelectrons per sensor system resolution 
element per sample time. 

(7) Apply a Poisson noise process to the array 
to simulate the noise effects of the 
photoelectron process. 

(8) Apply resolution degradation of the sensor 
and its associated electronics by Fourier 
transform techniques. 

(9) Add additional noise to simulate preamp and 
other system noise. 

(10) Discretize the image data to specific 
levels to simulate analog to digital 
conversion. This simulates the final 
recorded image. 

Figure 3 shows the results of a system simulation process 
for the same Images as shown in Fig. 2. 

Discussion 

The ability to simulate images of three-dimensional objects 
has been found to be a very useful tool for providing precisely 
defined and precisely variable images for use in system and 
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mensuration studies. Another application of object simulation 
which has yet to be fully explored is its use in image 
interpretation. (Ref. 4) Given the degraded image of a real 
object, the object and its image would be simulated. The 
simulation would be based on all previous knowledge of the 
object, the known lighting and viewing conditions, the 
characteristics of the imaging system, and information extracted 
from the actual degraded image. A comparison would then be made 
between the actual and simulated degraded image. Differences 
between these images which exceed the noise level of the real 
image would then indicate changes to be made to the assumed 
object. An iterative process would allow the true object to be 
constructed within the limits imposed by resolution and noise in 
the real image. 



M*** v 


- 13 - 




Node list: (x^Zj), . . . .(XgygZg) 

Connectivity list: (1,4,3,2), (4,8,7,3), (8,5,6,7), etc. 


(a). Representation of an object by polygons. 



FIGURE 1 
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§ Rome Air Development Center 2 

n MDC plans and executes research, development, test and 2 

• selected acquisition programs in support oi Command, Control Q. 

ft Communications and Intelligence (C 3 I) activities. Technical ? 
\ rnd engineering support loithin areas oi technical competence $ 
ft ** provided to BSD Program OUices (POs) and other BSD $ 

D elements. The principal technical mission areas are \ 

k communications, electromagnetic guidance and control, sun- & 
J veitlance o i ground and aexoipace objects, intelligence data jJ 
% collection and handling, information system technology, w 

* ionospheric propagation, solid state sciences, microtoave 5> 

i physics and electronic reliability, maintainability and & 
K compatibility. <ft 
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